The Balitsky-Kovchegov equation in full momentum space 
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We analyse the Balitsky-Kovchegov (BK) saturation equation in momentum space and solve it 
numerically. We confirm that, in the limit where the transverse momentum of the incident particle 
k is much bigger than the momentum transfer q, the equation admits traveling-wave solutions. We 
extract the q dependence of the saturation scale Q a (Y) and verify that Q S (Y — cste) scales as 
max(g,QT), where Qt is the scale caracterizing the target. 
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O ' I- INTRODUCTION 

o : 

The problem of understanding scattering near the unitarity limit in perturbative QCD has received a lot of attention 
over the past ten years. In the Regge limit of perturbative QCD, i.e. when the centre-of-mass energy in a collision 
is much bigger than the fixed hard scale of the problem, parton densities inside the projectiles grow with increasing 
energy, leading to the growth of the scattering amplitudes. As long as the densities are not to high, this growth 
is described by the BFKL equation [lj. But as the parton density becomes higher and the scattering amplitude 
] approaches the unitarity limit, one enters in a regime called saturation 0,0, where BFKL evolution breaks down. 
To study scattering amplitudes in the Regge limit, the gg-dipole model has been developed. This formalism 

constructs the light-cone wavefunction of a dipole (a quark- antiquark pair in the colour singlet state) in the leading 
00 • logarithmic approximation. The size of the dipole provides the hard scale that justifies the use of perturbation theory. 

As the energy increases, the dipole evolves and the wavefunction of this evolved dipole is described as a system of 
elementary dipoles. When this system of dipoles scatters on a target, the scattering amplitude has been shown to 
obey the BFKL equation. 

Interestingly enough, the dipole formalism was shown to be also well-suited to include density effects and non- 
linearities that lead to saturation and unitarization of the scattering amplitudes. Indeed, Kovchegov ja| derived an 
equation for the dipole scattering amplitude that reduces to the BFKL equation when the amplitude is small and that 
contains high-energy unitarity effects as the amplitude reaches unity. Although this equation is not exact, as it is a 
mean field approximation of the more complete infinite hierarchy of equations established by Balitsky , it has the 
advantage of being a closed equation for the scattering amplitude. This equation, which is usually referred to as the 
i-C BK (Balitsky-Kovchegov) equation, is propably the simplest and most accurate equation one can consider to describe 
the saturation regime of QCD. 

The BK equation is an evolution equation for A/"(r, b, Y), the imaginary part of the scattering amplitude of a 
dipole of transverse size r at impact parameter b and where the rapidity Y is the logarithm of the centre-of-mass 
energy. It is a non-linear equation that resums QCD fan diagrams in the leading-logarithmic approximation 
Many interesting features of this equation have been understood analytically 0, for the r and Y dependences. In 
particular, it was shown that the asymptotic solutions have the so-called geometric scaling |l(J property Af(r, h ~ 
0,Y) = Af(\r\QT^ s (Y)), where Qt is a scale characterising the target and fl s (Y) = exp(vY) is the dimensionlcss 
saturation scale. One has not been able however to extract any useful information on the b dependence. The analysis 
of the BK equation in impact-parameter space (see e.g. |lllll2j|) even shows a contradiction with confinement as the 
large-b dependence of the solutions develops a power-law tail with increasing rapidity. 

It was recently proposed [13j to analyse the equation in momentum space, by Fourier transforming A/"(r, b, Y), into 
A/"(k, q, Y) where k can be interpreted as the incoming transverse momentum of the dipole and q as the momentum 
transfer. It was shown in that the three variables k, q, and Y give information on the approach to the unitarity 
limit. In particular, it was predicted that the geometric scaling property should also be present in momentum 
space when |k| ^S> |q| as one should have asymptotically J\f(h,q,Y) = Af[\k\/ (\q\£l s (Y))]. The goal of this paper is 
to formulate and analyse the BK equation in momentum space, to solve it numerically and compare our results 
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FIG. 1: Linear contributions to the evolution of the dipole density in momentum space: (a) real gluon emission, (b) virtual 
gluon emission. 



with those predicted. The equation in momentum space has already been studied numerically in the b-independcnt 
situation |14j . Our purpose here is to extend those analysis to full momentum space including the transfer momentum 
dependence. 

The plan of the paper is as follows. In section [HJ we write the BK equation in momentum space and derive some 
analytical properties. In section ITTT1 we briefly explain the properties of asymptotic solutions and how they translate 
into geometric scaling. In Section llVl we give details about the numerical techniques used to solve the equation and 
present our results. Section W\ concludes. 



II. THE BK EQUATION IN MOMENTUM SPACE 

The BK equation in coordinate space reads 



2tt 



d 2 z 



(x-y) 2 



(x-z) 2 (z-y) 1 



[AA(x, z) + AA(z, y) - Af(x, y) - Af(x, z)M(z, y)] 



(1) 



where a = a s N c /Tr. The two-dimensional vectors x and y are the transverse positions of the quark and antiquark 
composing the dipole of transverse size r = x y at impact parameter b = (x + y)/2. The rapidity Y is not men- 
tioned explicitly in the amplitude but AA(x, y) always carries a Y dependence. Let us define the following Fourier 
transformation: 



#(k,q) = 7^1 fxfye^i-V*^^ 
(2tt) 2 J (x-y) 2 



(2) 



Then the equation verified by Af(k, q) can be written as 



dyrfik, q) 



d 2 k' 



(k-k') 2 



A/-(k',q)- 



,2 1 



(q-k) 2 _ 

(q-k') 2 T k' 2 



A^(k,q) 



a 
2^ 



d 2 k' Af(k, k')AA(k-k', q-k'). (3) 



This equation involves three contributions. The first integration contains the linear terms corresponding to real and 
virtual gluon emissions (see Fig^). The last term is the nonlinear correction which, as represented in Fig[3 arises 
from the resummation of fan diagrams. A trivial change of variable in the virtual correction term in Eq. © can put 
the evolution equation under a form emphasizing more clearly the cancellation of divergences between the real and 
virtual gluon emissions. We get 



dyA7(k,q) = - 



d 2 k' 



7T J (k-k 



A2 



A^(k',q)-- 



(q-k) 2 



k' 2 + (k-k') 2 (q - k') 2 + (k - k') 



Af(k,q) 



^ / d 2 k'rf(k,k')Ar{k-k',q-k'). 



(4) 
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FIG. 2: Non-linear contribution to the evolution of the dipole density in momentum space. This correspond to the resummation 
of fan diagrams. 



Let us now look at the properties of this equation in some particular cases. If the amplitude (in coordinate space) 
only depends on the size of the dipole i.e. 7V(x, y) =JV(x— y), then M(k, q) = 2n5^ (q)7v (k), where 



JV(k) 



d x 
2vrx 2 



W(x) 



obeys 



9yAA(k) = - 



d 2 fc' 



(k-k') 



>Y2 



AA(k') 



k' 2 + (k - k') 2 



,AT 2 (k). 



(5) 



(6) 



One recovers the usual equation discussed for b-independent situations. 

Let us now look at the equation in the forward limit q — * 0. Due to the fact that the non-linear term in is off- 
diagonal, the equation for the forward dipole amplitude involves non-forward contributions. However, if we assume 
that the amplitude is peaked around |q| = and has the factorised form 



A/-(k,q)=<K|q|)^o(k), 



(7) 



the non-linear term can be simplified 



d 2 fc'AA(k,k')AA(k-k',q-k')«2 7 rAA 2 (k) / k'dk' cj) 2 (k') 



If we normalize the q-dependent factor in Q in such a way that 



k'dk' 2 (fc'), 



we obtain a closed equation for the rapidity evolution of A/o(k) of the form 



d 2 k' 



7T J (k-k') 5 



Tvo(k') 



(k-k 



iAft(k), 



which is equivalent to the equation © obtained in the b-independent situation. It is however important to realize 
that, while the simplification is exact for b-independent amplitude, it requires assumptions on the small-q dependence 
of the amplitude in the forward limit. 



III. TRAVELING WAVES AND GEOMETRIC SCALING 



In this section, we briefly summarize the properties of the asymptotic solutions of the BK equation and replace this 
in the context of the new form J2J . 
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Let us start with the case of the b-independent equation ©. It has been shown p] that, in the diffusive approxima- 
tion, the equation for Af(k = |k|) is equivalent to the Fisher-Kolmogorov-Petrovsky-Piscounov (F-KPP) equation [l5| . 
This equation, well known in statistical physics, admits traveling-wave solutions at large rapidities and momentum 
k, which means that the solution is a function of Y — u c log(fc) only, where v c is the speed of the wavefront. More 
generally, if the evolution equation satisfies the following conditions: 

• M = is an unstable fix point and the equation contains a linear term, leading to a growth of the amplitude, 
and a non-linear damping term; this is realized in Q and J3J , 

• the initial condition is steep enough, which in the QCD case is ensured by colour transparency, 

• the linearised equation admits a superposition of waves as solution 

rc+ioo , 

M Vm (k,Y) = -^a ( 7 )exp[c(7)r-7log(fc 2 /fc 2 )] , (8) 



then the asymptotic solution of the full equation takes the form 



Af(k,Y) 



k 2 



1 ~Tc 



(9) 



for k 3> fco where the saturation scale Q S {Y) grows as ko e Vc . The critical exponent 7 C and the critical speed v c are 
given by 



w'(7), 



meaning that the group velocity is equal to the phase velocity. 

In the case of the BK equation, the corresponding linear equation is the BFKL equation, whose solutions have 
been extensively studied jl6L Il7l Il8l Il9j . For the impact-parameter- independent case, Munier and Peschanski have 
proven that the three conditions for traveling waves are fulfilled. In a recent study, we have shown 0] that this 
can be extended to the full BK equation including momentum transfer. In this case, the reference scale k becomes 
proportional to the typical scale of the target Qt or to the momentum transfer |q| whether |q| is smaller or larger 
than Qt- Before going into the numerical studies, we recall the important points concerning geometric scaling at 
non-zero momentum transfer. 

In the previous work one of the key points was to show that the BFKL solutions can be put under the form © 
when |k|^> |q|. Since the definition J2J of the Fourier-transformed amplitude introduces an additional (x— y) 2 w.r.t. 



the usual one, the solution is slightly different to the one we used in ^ 
we shall see now. Let us start with the well-known result [l6l IT^. ITsl u\ 



. However our conclusions are unchanged as 
that 



Min(x,y) 



dry 
2in 



o a x(l)Y 



x y 



2l 



00 (7) 



is solution of the BFKL equation in coordinate space with 



x(7)=2V(l)-^(7)-^(l-7)- 



(10) 



and where </>o(7) is a function that specifies the initial condition. Note that we kept only the dominant component at 
high energy which corresponds to vanishing conformal spin n = 0. If we now Fourier transform this solution according 
to one obtains 



Min(k,q) 



_L e a ^r(k, q )<M 7 ,q) 



with 



/ 7 (k,q) 



r 2 (7) 



r 2 (i + 7 ) \k\ 



Q_ 

4k 



27-1 



2F1 (7, 7; 27; |j 2 Fi ( 



7, 7; 27; r 



(7-*- 1- 7). 



(11) 



A factor antisymmetric under the replacement 7^1 — 7 and as been absorbed in </>o(7, q). Note that in (|ll|l we use 
the complex representation of vectors: x — xi+ix2 and x = x\ — 1x2 for the vector x = (xx,X2); however from now 
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on, we shall use x to denote the modulus of the vector. Although this is not straightforward, one can explicitly check 
(see Appendix [SJ that is an eigenfunction of the linear kernel of eq. with the correct eigenvalue (|10|) . When 
k q, the hypergeometric functions go to 1 and / 7 (k, q) becomes a simple power of q/k. The full solution of the 
linear part of © thus becomes 



Mi„(k,q) = 



dj 

2iTT 



Ml, q) exp [ax{l)Y - -f\og{k 2 / q 2 )} , 



where irrelevant factors have again been absorbed in the initial condition 0o(7,q)- This expression proves that the 
conditions for obtaining traveling waves and geometric scaling arc fulfilled. More precisely, including all the terms 
common to each evolution equation of this universality class |9j, it predicts 



Q 2 S (Y) = q 2 n 2 s (Y) 



q 2 exp 



7c 2 7c 7^ 



2tt 



1 



c*X" (7c) VY 



0(1/Y) 



together with the leading edge of the traveling front 



■Af(k,q) 



Tc (q) log 



q 2 n 2 (Y) 



q 2 n 2 (Y) 



exp 



1 



2a X "(7c)^ 



log 2 



q 2 n 2 (Y) 



(12) 



(13) 



When q becomes smaller that the typical size Qt of the target, one can show that Qt substitutes to q in the saturation 
scale. Q s is thus expected to be proportional to q only at large q. The value of 7c is 0.6275. Note that those behaviours 
are the same as the ones predicted by the forward analysis. 

Let us also point out that all these results have been studied in terms of the momentum transfer q and could not 
be derived in terms of the impact parameter b. One way to understand this is to look at the BFKL equation. One 
clearly sees that the equation in momentum space (see (|2Jl) is local in q while the equation in coordinate space (see 
couples different values of the impact parameter b. In the following section, we shall explicitly check that these 
predictions are observed by numerical simulations of the full equation . 



IV. NUMERICAL ANALYSIS 



A. Description of the method 



Instead of working directly with the momentum k, we shall use the variable p 
conjugate of r. The Fourier transform of the dipole amplitude then becomes 



k — q/2, which is the canonical 



■#(p;q) 



1 



(27T) 5 



d 2 r 



d 2 be lr p e lb ci Af(r- b), 



and the BK equation in momentum space takes the form 

d 2 p' 



cVAA(p;q) = 



2nJ (p'-p) 2 



2AA(p';q)- 



(P + f) 



(p-ir 



p'+§) +(p'-p) 2 (p'-i) +(p'-p) 2 



AT(p;q) 



2tt 



— / d y 'N p - — + -; - + p AM P - — - T ; - - p 



4' 2 



q q 



4' 2 



(14) 



Although this formulation seems more complicated, it emphasises more clearly the symmetry properties of the BK 
equation under the transformation k > q — k which simply becomes p — * — p in the new variables. Let us notice 
that the solution of the linear problem can be seen as depending only on the ratio p/q. Also, the limit fc 3> q we are 
interested in is equivalent to p ^> q. 

For numerical studies, we shall therefore use L = log(Rp), L q = log(Rq), where R is the typical target size, and (f>, 
the angle between p and q. It is obvious from 1|14|) that we have 

Sf(L, L q , 0) = N{L, L q , -0) - Af(L, L q , tt-4>)= Af(L, L q , n + 0), 



provided the same equality holds for the initial condition. In practice, we have computed the amplitude M{L, L q , (f) 
on a lattice of n p x n q x points regularly spaced with L min < L < L max , L q min < L q < L q maK and < 4> < \. 
We have then discretised the integral in the r.h.s. of i|14[l and solved the resulting set of differential equations by a 
Runge-Kutta method. 
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B. Results 

In this section, we shall study the properties of the numerical solutions of equation lfH|) . We have studied two 
choices of initial distribution: 



1 

Y=0 p 2 + R~ 



#(p;q)„ „ = 2 , D -2 ( 15 ) 



and 

A^( P ;q) 



Y=o p 2 + R- 2 1 + q 2 R 2 ' 



1 (16) 



Although the first choice seems more natural, the exponential damping at large q leads to very small amplitudes 
producing numerical results difficult to handle. Therefore, in order to make the study at large q easier, we shall also 
consider the less-steep behaviour The value of the coupling a has been fixed to 0.2. We shall not consider the 

effect of the running of the strong coupling constant in this paper. Before presenting the results, let us notice that 
the dependence of AT on the angle tfi is not really interesting for our purposes, hence, we shall concentrate on the 
amplitude averaged over <f>. Concerning the choice of the discretisation parameters, we have used 101 values of L 
between —5 and 20, 53 values of L q between —5 and 8 and 9 values of <j> between and 7r/2. 

To begin with, we shall check if the traveling waves predicted for k 3> q in are indeed observed in numerical 
simulations. This is done by plotting the rapidity evolution of the dipole amplitude as a function of k (or, equivalently, 
p) for different values of q. As clearly seen in Figure |3 as the rapidity increases, a traveling wave is formed at large 
values of p. For these results, the initial condition used is ((TBI) . 

To observe this more precisely, it is interesting to study the properties of the saturation scale Q s . In it was 
predicted that Q s should evolve asymptotically like a ^-independent power of the energy, and, at fixed rapidity, should 
grow linearly with q when q becomes larger than the inverse size of the target. Let us first check the high-energy 
dependence of the saturation scale. Due to the g-dependence of the initial condition (|16(l . one cannot simply extract 
the saturation scale by solving N(p = Q s (Y,q);q) — Nq since, for sufficiently large values of q, the amplitude will 
always be smaller than Nq. To fix this problem, we rely on the fact that the asymptotic behaviour of log(Q s ) is 
expected to be insensitive to the choice of Nq. We shall thus use a ^-dependent threshold No(q), taking into account 
the damping at large q introduced by the initial condition. This leads us to solve 

JV(p = Q s (Y,q):,q) = #„(?) = jf^' ( 17 ) 

for different values of Y and q. In Figure Ufa), where we have chosen iVo(O) = 0.05, we show that whatever q is, 
log(Qs(Y)/Q s (Yo)), with Yq fixed to 1, goes to a constant value at asymptotic rapidities. In fact, one can even check 
that our numerical results are consistent with the expansion 

d Y \og [Q 2 s (Y,q)] = 

valid at large Y. Notice that, in this expression, the g-dependent non-universal term C2Y~ 2 is required in the rapidity 
range under consideration. 

The extraction of the g-dependence of the saturation scale is more tricky since the definition (|17(l is obviously not 
sufficient. Basically, one has to obtain the position of the wavefront in a way that does not depend on any threshold Nq 
at all. Our method consists in studying the so-called reduced front p 2jc J\f(p, q) which is expected to show a maximum 
in p in the region where geometric scaling holds [jj- However, since the violations of geometric scaling in the ultra- 
violet appear outside of our lattice, we shall define Q s as the maximum of a modified version of the reduced front 
pAf(p,q) i.e. replace 7 C = 0.6275 by 0.5. Physically, the maximum of this quantity can be considered as a measure 
of the momentum for which the saturated regime is replaced by the traveling front. We have applied this method to 
extract the average 1 saturation scale for 20 < Y < 25 and checked that the result does not change if one replace 0.5 
by another value close to it. The result is shown in Figure Ufb) and we see that the saturation scale obtained from 
this method is constant for small momentum transfer and shows a linear increase as q becomes larger than R . For 




1 We average over an interval of rapidity in order to minimise numerical uncertainties. 
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FIG. 3: Rapidity evolution of the dipole density as a function of p for different values of q. For each plot, we show the amplitude 
for Y varying between and 25 by steps of 2.5. 
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FIG. 4: This figure presents the results concerning the saturation scale: (a) shows that the evolution of the saturation scale 
with Y has the predicted behaviour and (b) represents the g-evolution of the saturation scale (see text for details). 
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FIG. 5: This figure shows the Y dependence of the amplitude 7V(p;q). In the left plot we have fixed log(i?q) = —3, while 
log(i?g) = in the centre plot and log(Tig) = 3 in the right one. 
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FIG. 6: This figure shows the Y dependence of the amplitude A/(p; q) compared to the amplitude obtained by BFKL evolution. 
In each pair of curves, the upper one is obtained from BFKL evolution and the lower one includes nonlinear effects. In the left 
plot we have fixed \og(Rq) = —3, while \og(Rq) = in the centre plot and log(Rq) — 3 in the right one. 
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FIG. 7: Rapidity evolution of the dipole density as a function of q for different values of p. In the left (resp. centre, right) plot 
we have set log(Rp) = —3 (resp. 0, 3). 

very high values of q (\og(Rq) > 2), the results start to diverge from this proportionality. We attribute this feature 
to the fact that the asymptotic traveling-wave behaviour is not yet reached, and at those high values of q, we are not 
extracting the true saturation scale. 

Let us say a few words about the evolution in the regime fc~g. As shown in Figure^ this region presents oscillations 
and the amplitude seems to saturate to a constant value at large Y. To see more precisely the effect of the nonlinear 
term in we can compare the y-evolution of the amplitude presented in figure with the results obtained from 
BFKL evolution i. e. keeping only the linear term in @ . The comparison of these two evolutions is presented in Figure 
|SJ This figure shows interesting features. First, one clearly sees the damping due to the inclusion of the nonlinear 
effects. Then, one observes the exponential growth of the BFKL amplitude with rapidity. Finally, we see that at high 
energy and for a given value of q, the BFKL amplitude becomes degenerate for k ^ q. This feature, closely related 
to the plateau observed in the infrared region in Figure [21 shows that q acts like an infrared cut-off p0| . 

Finally, let us see how the q behaviour of the amplitude evolves with rapidity. In Figure \7\ we show how the 
initial condition (|15fl evolves with rapidity for different values of p, as a function of q. First of all, it appears that 
the extension in q remains roughly constant as the rapidity increases. Then, we see an increase of the amplitude for 
small values of q. As noticed in [13J, the fact that the q-behaviour of the amplitude does not develop a power-law tail 
when the energy increases as it is the case with impact parameter, is one of the main advantage in studying the BK 
equation in momentum space. 

V. CONCLUSION 

Let us summarize the main results of this paper. First, we have proposed a Fourier transform of the dipole scattering 
amplitude which leads to a new form of the BK high-energy evolution equation. This new form involves the amplitude 
in momentum space instead of coordinate space (see equation ©). We have shown that the momentum-transfer- 
dependent equation reduces to the usual momentum-space evolution in the ^-independent case. This equation is also 
obtained in the forward case q = 0, assuming a steep q dependence of the amplitude. 

We have then checked numerically the properties of the asymptotic solutions of the proposed q-dependent BK 
equation. We have shown that all the properties predicted theoretically in jl3| are indeed reproduced by numerical 
simulations. More precisely, one obtains traveling waves at large rapidities, leading to geometric scaling in the region 
k ^> q. The speed of the traveling front converges to the g-independent value predicted by previous analysis. Concern- 
ing the saturation scale, we have shown that it goes from a constant behaviour to a linear increase with momentum 
transfer when q goes from values smaller than the target typical scale to larger values. While the coordinate-space 
version of the BK equation leads to a diffusion in impact parameter already at the level of the linear kernel, the 
corresponding equation in momentum space remains local in q. This means that, to the contrary of what happens in 
impact-parameter space where the evolution generates a perturbative tail at large &, the q behaviour does not show 
an increase at large values of q. This important property is also the reason for which momentum space is the natural 
frame to discuss traveling-wave solutions and geometric scaling. 

In addition to these results, the possibility to solve the full BK equation allows for various extensions. From 
the phenomenological point of view, it can help to test saturation effects in many processes including e.g. diffractive 
photoproduction of vector mesons. On more theoretical grounds, one may be tempted to analyse the odderon equations 
Ref. [21| . Due to the parity of the odderon, this contribution cannot be introduced at the level of the b- independent 
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BK equation. In all these cases, it would be interesting to introduce the effect of the running coupling. Finally, this 
equation can also serve as a basis for studies of fluctuation effects and the corresponding Langevin equation derived 
very recently |2^ . 
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APPENDIX A: SOLUTIONS OF THE LINEARISED EQUATION 

In this appendix, we write down the main lines that show that (|llf> is a solution of the linear part of equation (|3J) 
with the correct eigenvalue. We start by going to the complex representation and use the integral representation of 
/ 7 (k, q), see Appendix B of 0. One has, up to a 7-dependent factor: 



p{ x ) = \ x \ 2 ~< I d z v M Z7_Z |1 - \l-xv\ 



27-2 



127-2 



-27 



with x = q/k. If we perform the change of variables k' — ► v = k! /k in (the linear part of) 0, we are left with the 
integral 



Z7T 



d 2 V 



\l-v\ 



2P(x/v) - 



a 

+ — 

47T 



l-x 



\i — v\ y \v\ \v — x\ £ 



The second integral can be computed, e.g. by dimensional regularisation. We also perform the change of variables 
v — * u — 1/v in the first integral, leading to 



P{x) 



d 2 U 

\l-u\ 



1 



r(ux) 



1+ 1-tt 



log 



1 



Inserting the integral representation of / 7 (a;), one obtains 

d 2 u 



P(x) 



-\x^ 

7T 



\x\ zl I d 2 v \v\ ZJ - z \1 - v\ 2 *>- 2 lim 



6^0 



II - u 



2-2s 



| U | 2 T 2 |1 — VXU\ 27 



d 2 u 



|1 — vx\ 



-27 



\l-u\ 2 - 2e l + |l-d 2 



a 1 

— log 

2 B 



1 



l-x' 



where we have explicitly regularised the divergence at u = 1. 

Performing the integration over u and taking the limit e — > 0, we get 

P{x) =a X { 1 )P{x)+aW{x), 



with 



•log 



1-x- 



+ |.t| 27 / d 2 v | V | 27 - 2 |1- V | 27 - 2 [(l~xvyW £2 F 1 ( llT , 1 + e;vx) + (l-xv)^d e2 F 1 ( 1 , Tll + e;vx)] e=Q . 



The derivative of the hypergeometric function can be explicitly computed and we are finally left with the integration 



i? 7 (z) 



1 27 



,2 1 127—2 1 1 |27-2 M _ 1-27, (\1 — XV\ 

a v \v\ 1 — v\ 1 — xv\ log 1 



\l-x\- 
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Finally, the conformal transformation v — frEffr gives 

R^x) = ^f - [d 2 u |u| 2 ^ 2 |1 - u\^- 2 |1 - cnx^ log (ji^C) = 
* J \^ 1 1 — xu\ ) 

Therefore, hV(x) = and f f (x) is an eigenfunction of the linear kernel with xil) as associated eigenvalue. 
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